Multi-Species Asymmetric Exclusion Process in Ordered 

Sequential Update 



M.E.Fouladvand" 6 Q , F. Jafarpour a 6 Q 



a Department of Physics, Sharif University of Technology, 
P.O.Box 11365-9161, Tehran, Iran 
b Institute for Studies in Theoretical Physics and Mathematics, 
P.O.Box 19395-5531, Tehran, Iran 

Abstract 
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Product Ansatz (MPA), we obtain the relevant algebra, and study the uncorrelated 
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1 Introduction 



One dimensional models of particles hopping in a preferred direction provide simple non- 
trivial realizations of systems out of thermal equilibrium [1,2,3,4]. In the past few years 
these systems have been extensively studied and now there is a relatively rich amount of 
results, both analytical and numerical, in the literature, (see [1,4] and references therein). 
These types of models which are examples of driven diffusive systems, exhibit interest- 
ing cooperative phenomena such as boundary-induced phase transition [5], spontaneous 
symmetry breaking [6,7] and single-defect induced phase transitions [8,9,10,11,12] which 
are absent in one dimensional equilibrium systems. 

A rather simple model which captures most of the mentioned features is the Asymmetric 
Simple Exclusion Process (ASEP) for which many analytical results have been obtained 
in one dimension [1,4,13]. Besides its usefulness in describing various problems such as 
kinetic of biopolymerization , surface growth , Burgers equation and many others (see 
[4] and references therein), ASEP has a natural interpretation as a prototype model de- 
scribing traffic flow on a one-lane road and constitutes the basis for more sophisticated 
traffic flow models [14,15,16]. 

Derrida et al were first to apply Matrix Product Ansatz (MPA) in ASEP with open 
boundaries [17]. Since then, MPA has been applied to many other interesting stochastic 
models such as ASEP with a defect in the form of an additional particle with a different 
hopping rate [11], the two species ASEP with oppositely charged particles moving in 
same (opposite) directions [6,12,18] and many others. MPA has also been shown to be 
successful in describing disorderd ASEP- like models. Evans [19] considered a model on a 
ring where each particle hops with its own specific rate to its right empty site if it is empty 
and stops otherwise. This model shows two phases. In low densities the hopping rate of 
slowest particle determines the average velocities of particles (phase I). When the den- 
sity of particles exceeds a critical value, it is then the total density which determines the 
average velocity and the slowest particle looses its predominant role (phase II). In spite 
of its many nice features both theoretically - drawing a relevance between Bose-Einstein 
condensation and the observed phase transitions - and idealistic - a better modeling of a 
one way traffic flow - the possibility of exchanging between particles overtaking has not 
been considered. This possibility is crucial in describing a more realistic traffic flow model. 

Very recently in [20], a multi-species generalization of ASEP has been proposed such 
that exchange processes among different species has already been implemented. In this 
model, there are p-species of particles present in an open chain with injection(extraction) 
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of each species at boundaries. Each particle of i-type ( 1 < i < p ) hops forward with 
rate Vi and can exchange its position with its right neighbour particle of j-type with rate 

Vi - Vj. 

The subtractive form of exchange rates allows that only fast particles exchange their 
positions with slow ones. As can be seen, this is a more realistic model for traffic flow 
in which fast cars pass slow ones in a two-lane one-way road. In [20] using MPA, an 
infinite dimensional representation of the quadratic algebra is obtained but the form of 
currents and density profiles could not been obtained by this infinite dimensional rep- 
resentation. Instead, the simple case of one dimensional representation was considered. 
Although restricting the algebra to be one dimensional, will cause to loose all the corre- 
lations, but still many interesting features such as a kind of Bose-Einstein condensation 
and boundary induced negative current [21], appear even in this simple uncorrelated case. 

Most of the above mentioned models have been defined in continuous time, where the 
master equation of the stochastic process can be written as a schrodinger- like equation 
for a " Hamiltonian" between nearest-neighbours [4,22]. In contrast, one can use discrete- 
time formulation of such random processes and adopts other type of updating schemes 
such as parallel, sub-parallel, forward and backward ordered sequential and particle or- 
dered sequential (see [23] for a review). The MPA technique has been extended to a 
sublattice parallel updating scheme [24,25,26] and in the case of open boundary condi- 
tions, to ordered sequential scheme [27,28]. Although in traffic flow problems, parallel 
updating is the most suitable one, only few exact results are known [15,29,30]. 

In general, it's of prime interest to determine whether distinct updating schemes 
can produce different types of behaviour. The present analytical results show that with 
changing the updating scheme of the model, general features and phase structure re- 
mains the same but the value of critical parameters may undergo some changes. In [23], 
Schreckenberg et al have considered ASEP under three basic updating procedures. Sim- 
ilarities and differences have fully been discussed. Evans [29] has obtained analytical 
results in ordered and parallel updates for his model which was first solved in random 
sequential updating in [19]. He has demonstrated that the phase transition observed in 
[19] persists under parallel and ordered sequential updating. 

In this paper, we aim to study the p-species model introduced in [20] under ordered 
sequential update scheme and will show that the features observed in [20] are reproduced 
in ordered updating as well. Our results will be reduced to those of [23] when we set p=l. 
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The organization of the paper is as follows. In section 2, we briefly explain the p- 
species ASEP with random sequential updating and then describe the MPA in backward 
ordered sequential updating for jo-species ASEP and will obtain the related quadratic 
algebra. Section 3 contains the mapping of algebra of section 2 to that of [20] and 
includes the expressions for the currents and densities of each type of particles in the MPA 
approach. Section 4 is devoted to the one dimensional representation of the quadratic 
algebra and the infinite limit (p — > oo) of the number of species. In this limit, we use a 
continum description of current-density diagrams of the model. In section 5, we consider 
the forward ordered updating and discuss the similarities and differences between forward 
and backward updating. In contrast to the usual ASEP where particle-hole symmetry 
allows for a map of result between forward and backward updating [23], here we don't 
have particle-hole symmetry and hence should separately consider the forward updating. 
At the end of this section, we discuss the intimate relationship between sub-parallel 
scheme and ordered sequential [32]. Section 6 concludes the model with ordered updating 
on a closed ring. We obtain current-density diagrams for both backward and forward 
updating. The paper ends in section 7 with some concluding remarks. 

2 The Model 

2.1 p-species ASEP in ordered sequential updating 

In this section we first briefly describe the p-species ASEP introduced in [20]. This 
model consists of a one dimensional open chain of length L. There are p species of particles 
and each site contains one particle at most. The dynamics of the model is exclusive and 
totally asymmetric to right. Particles jump to their rightmost site provided that site is 
empty, time is continuous and hopping of a particles of type i (1 < i < p) occurs with 
the rate V{. To cast a more realastic model for describing traffic flow, there has been 
considered the possibility of exchanging of two adjacent particles i.e. two neighbouring 
particles of types (j) and (i) swap their positions with rate Vj — Vi , Vj > Vi. This 
automatically forbids the exchange between low-speed and high-speed particles so it's 
a natural model for a one way traffic flow where fast cars can overtake the slow ones. 
Denoting an z-type particle by Aj, and a vacancy by 0, the bulk of the process is defined 
by: 

Ai<f> — > with rate (i — l,...,p) (1) 

AjAi — > AiAj with rate Vj — Vi {j > i — 1, ...,p) (2) 
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In order for all the rates to be positive, the range of Uj's should be restricted as: 



vi <v 2 < v 3 .... < Vp (3) 

To complete the process, one should consider the possibility of injection and extraction 
of particles at left and right boundaries. The injection (extraction) of particles of type i 
at left (right) boundary occurs with the rate (Pi). 

This completes the definition of the model. Denoting the probability that at time t, the 
system contains particles of type r« (ti = refers to vacancy) at site % (0 < r« < p, 1 < 
i < L) by P(ti, r 2 , r^, t), one can write the stationary state P s (ti, T2, Tl) in form of 
a Matrix-Product-State (MPS) 



P s (n,...,T L )~ < W\D T1 ...D TL \V > (4) 

in which D n (0 < Tj < p) is an ordinary matrix to be satisfied in some quadratic alge- 
bra induced by the dynamical rules of the model and the vectors \V > <W\ (reflecting 
the effect of the boundaries) act in some auxiliary space [31,32]. Denoting D by E, the 
quadratic algebra reads [20] 



DiE = -A + E (l<i<p) (5) 

DjDi = ]_ A viDj - VjDi) (l<i<j<p) (6) 

The vectors \V > and < W | satisfy 

Di \V>= V j\V> (7) 

<W\E =<W\— (8) 
pcti 

The rest of [20] concerns with obtaining densities and currents of species from the 
above relations. We will come back to these results while comparing them with ours 
in the comming sections. In what follows, we describe p-species model under ordered 
sequential update. 

As stated in the introduction, in ordered sequential updating, time is discrete and the 
following events can happen in each time-step 



Ai4> — > <f)Ai with probability (i — l,...,p) (9) 
AjAi — > A { Aj with probability fa (j > i — 1, (10) 



We do not fix the form of fjiS and as will be seen, they will be fixed later. Particles are 
also injected (extracted) at the first (last) site with the probability a.i (Pi). We denote the 
probability of the configuration (t±, tl) at N'th time-step by P(ti, ...,tl',N). We make 
a Hilbert space for each site of the lattice consisting of basis vectors {|r >, r = 0, ...,p} 
where |r > denotes that the site contains a particle of type r (vacancy is a particle of 
type 0). The total Hilbert space of the chain is the tensor product of these local Hilbert 
spaces. With these constructions, the state of the system at the N'th time-step is defined 
to be \P, N > so that 

P(t u ...,t l ;N) :=< Ti ,...,t l \P;N> (11) 

In ordered sequential updating one can update the system from right to left or from 
left to right. In general these two schemes do not produce identical results, so it is 
necessary to consider both of them separately. We first consider updating from right to 
left (backward). The state of the system at (j + l)'th time-step is obtained from j'th 
time-step as follows 

\P,j + l>=T^\P,j> (12) 

where T<_ is 

T<- = £1^1,2^2,3 ■ ■ ■ Tl-2,l-iTl-i,lRl (13) 

with 

Li = L<g>l<g>---<g>l , R L = 1 <g> 1 <g> ... <g> R (14) 
T i)i+ i = l(8)l---J > (8)T(8)J > ---l(8)l (15) 

i-1 i+2 

According to (13), updating the state of the system in the next time-step consists 
of the L + 1 sub-steps. First the site L is updated: if it's empty it's left unchanged, 
but if it contains a j-type particle (1 < j < p) , this particle will be removed with the 
probability fa from the site L of the chain, then the sites L and L — 1 are updated by 
acting T l _x,l on |r i _ 1 > ®R\r L >. The effect of Tl-\,l is to update the site L — 1 and 
L according to the stochastic rules (9) and (10). After updating all the links from right 
to left, one finally updates the first site: if it's occupied it's left unchanged, if it's empty 
then a particle of type i (1 < i < p) is injected with the probability c^. This procedure 
defines one updating time-step. After many steps, one expects the system to reach its 
stationary state \P S > which must not change under the action of TU and therefore is an 
eigenvector of with eigenvalue one 
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\P S >=T^\P S > (16) 

The explicit form of T, R and L can be written as 

v v 
T = J2vi(E 0i ®E i0 - Eu(g)E 00 )+ ]T f ji (E ij ®E ji -E jj ®E ii ) + I (17) 

i=l j>i=l 

R = ir[3 l (E 0l -E ll )+I (18) 
i=i 

L = j^a t {E lf) -E m ) + I (19) 
i=i 

Here the matrices E^ act on the Hilbert space of one site and have the standard 
definition (E^m = 5 ik 8ji. 

2.2 Matrix Product Ansatz (MPA) for ordered sequential scheme (back- 
ward) 

In this section we introduce MPA for the p-species model with right to left ordered 
sequential updating scheme. As shown by Krebs and Sandow [31], the stationary state 
of an one dimensional stochastic process with arbitrary nearest-neighbour interactions 
and random sequential update can always be written as matrix product state (MPS) 
[31]. In [32] Rajewsky and Schreckenberg have genaralized this to ordered sequential 
and sub-parallel updating schemes which are intimately related to each other. Following 
[17,23] we demand that 

P a (Ti,...,T L )~ <W\D n ...D TL \V> (0<r t <p) 

where the matrices D ,...,D P and the vectors \V > , < W\ are to be determined. 
Let's first write the above MPS in a more compact form via introducing two column 
matrices A and A 
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(elements of A and A are usual matrices) so we formally write 

\P S >= — « W\A® A® ... <g> A\V » (20) 

where the normalization constant Z L is equal to < W\C L \V > with C = E + Yn=i Di. 
The bracket << • • • >> indicates that the scalar product is taken in each entry of the 
vector A® A ■ ■ ■ <g) A. One can easily check that (20) is indeed stationary i.e. T^\P S >= 
\P S > , if the following conditions hold 



RA\V >= A\V >, (21) 
T(A®A) = A®A, (22) 
<W\LA=<W\A (23) 

This simply means that a "defect" A is created in the beginning of an update at site 
j = L, which is then transfered through the chain until it reaches the left end where it 
disappears. Equations (17-19) and (21-23) lead to the following quadratic algebra in the 
bulk : 



EDi + ViDiE = EDi 



D l D j , 



, i = 


1,.. 


;P 


, i = 


1,.. 


;P 


, i = 


1,.. 


;P 


> i = 


1,.. 


;P 


> i = 


1,.. 


;P 



+ D i D j 
(1 - fji)Djbi 
and following relations 

v 

< W\(l-^2ou)E =< W\E 



(24) 
(25) 
(26) 
(27) 
(28) 



(29) 



i=l 
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< W\(aiE + Di) =< W\Di , i = l,...,p (30) 

(E + j2PiDi)\V>=E\V> (31) 
i=i 

(1 - Pi)Di\V >= D t \V> , i = l,...,p (32) 

3 Mapping of the p-species Ordered Sequential Al- 
gebra onto Random Sequential Algebra 

In this section we find a mapping between the algebra (24-32) and (5-8). This mapping 
for p = 1 (usual ASEP) was first done in [33] where it was shown that apart from some 
coefficians, ASEP in an open chain with either random or ordered update, leads to the 
same quadratic algebra. Here we show that this correspondence again holds for p-species 
ASEP. We first demand 



E = E + e (33) 
Di = Di- di , i = l,...,p (34) 

where e and di are c-numbers. Putting (33,34) into (24-32) one arrives at 



v i D i E = (l-v i )eD i + d i E , i = l,...,p (35) 
f ji D j D i = d j D i -d i (l-f ji )D j , j>i = l,...,p (36) 

< W\E =< W\e(- - 1) (37) 
a 

Di\V>=^\V> i = l,...,p (38) 

Pi 

in which a = Y?i=i a i an d the following constraints must be satisfied 



e = Vdj , ai = {-)di (i = l,...,p) (39) 
i=i e 

One should note that as soon as restricting the algebra (24-32) to the conditions 
(33,34), the probabilities of injection are no longer free and are restricted by (39). Up 
to now the exchange probabilities fji have been free, however we have not yet checked 
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associativity of the algebra (35,36). Demanding associativity fixes these exchange prob- 
abilities to be 

fa -7—— L > j>i = l,...,P (40) 

1 Vi 

Remark: according to the discrete-time nature of updating procedure, fj^s are more 
precisely, the conditional probabilities i.e. they express the probability of exchanging 
between j and z-type particles provided that the z-type particle does not hop forward 
during the sub time-step. Thus 

prob(- ■ ■ AiAj ■ ■ ■ ; N + 1| • • • A 3 A t ■ ■ ■ ; N) ~ /^(l - Vi ) = Vj - v t (41) 

Threfore we see that overtaking happens with a probability proportional to the the 
relative speed. With this requirement (35-38) yield 



DjDi 



Vj — V; 



ViDiE = (1 - Vi)eDi + diE , i = l,...,p 
{dj(l - vt)Di - di(l - Vj)Dj} , j>i = l,...,p 



< W\E =< W\e( 



1 



dj , 

A|K>=^>, 



a 

i = l,...,p 



(42) 
(43) 

(44) 
(45) 



(42-45) is the mapped algebra of p-species ASEP in backward ordered sequential 
updating onto random sequential updating. It can be easily verified that similar to one- 
species ASEP [17], any representation of the algebra are either one or infinite dimensional. 
In the following ZVs and E are explicitly represented 



E = 



/0 . . . 
10.. 
10. 
.010 
. . . 
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/ , XiO-Vj) x (l-^) 2 x (1-f,) 3 
Ai Ai — — — Ai — — — 

i 



(l-Vj) 1 (j-^if 1 







Vi 



V • 



where i] is a free parameter ( we have a class of representations). 

(46) 



with Aj = T , 

Using (45), we multiply both side of (43) on \ V > and we obtain 



Vj(l - Pi) - Vi(l - Pj) = f3j - A , j>i = l,...,p 
solving this equations yield 

A = (l+7K-7 , i = l,.-,p (47) 
in which 7 is a new parameter which can be written in terms of known quantities. 

TU{Pi-Vi) 



7 = 



Requiring that all the probabilities to be positive, leads to the following condition on 



Vi'S 



7 



7 + 1 



< vi < v 2 ... < v p < 1 , 7e[0, 00 [ 



(48) 



We conclude this section with formulas for the current operators. In contrast to 
random sequential updating where currents are local i.e. caused by at most a single 
hopping of particles, in the ordered sequential updating, the currents are highly nonlocal 
which to say can have many hoping sources according to the multiplicative nature of 
transition matrix T<_. In ordered sequential updating the mean current in the iV'th 
time-step through the site k is defined by 



< n k >N+1 - < n k >N~< Jk~l,k >N - < Jk,k+1 >N 



(49) 



Our attention is concentrated on the stationary state so N should go to 00. With 
introducing a bra vector 
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<S\-.= J2 <n,...,T L \ 

T1,...,TJV 

the l.h.s of (49) can be written as 

< 5|4 i) T < _T^|P(0) > - < S\n^T^\P(0) > (50) 
which in turn yields 

< nf > N+l - < nf > N =< S\[n®,T^]\P a > (51) 

We have used the fact that < S\T^ =< S\ which is justified if T<_ is the transfer matrix 
of a stochastic process . Evaluating the commutator in (51), everything is expressed in 
stationary state expectation values of densities which using MPS (20) would finally leads 
to the expression for the current of i-type particles from the site k — 1 to k 

w < w\c k -*j^c L - k \v , 

<^-i,*>-- <w \c L \v> {b ' 

in which 

JW = ViDtE + t V -^D l D 1 - ± ±—±1),!); (53) 

j>i 1 v 3 j>i 1 u » 

and 

C = E + j^D i (54) 
i=i 

The first term in (53) is due to hopping of the z-type particles, the second term cor- 
responds to the exchanges between an i-type and all the particles with lower hopping 
probabilities than it and finally the last term expresses the exchanging between all the 
particles with higher hopping probabilities and the z-type particle. 

Using (33), (34) and the bulk algebra (42) and (43) one easily concludes that 

J w = diC (55) 
So the current and density of z-type particles through (at) site k are respectively given 



by 

^ 7.w ^ =rl ;_ 

< w\c L \v > 
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Therefore all the currents are proportional to the average current , however 
has a nontrivial dependence on hopping probabilities. The next section is devoted to 
the one dimensional representation of the algebra (42-45). This case corresponds to the 
steady state characterized by a Bernouli measure. In spite of its simplicity, still some 
interesting features survive in a one dimensional representation. 



4 One Dimensional Representation and Infinite-Species 
Limit 

4.1 One dimensional representation 

The simplest representation of the algebra (42-45) is to take the dimension of the 
matrices to be one. For later convenience, let us replace all A's by ^ where p is 

the number of species. Denoting -y and E by c-numbers, ^ and S respectively, from 



equations (44) and (45) we have 



A = 7T-# ' * = (58) 

(1 + j)Vi - 7 a 

Putting these numbers in (42) leads to 

Vi = l or = 1 (59) 

a 7 

The case Vi — 1 corresponds to the ordinary 1-species ASEP which has been exten- 
sively studied. Using (47) the second condition can be written as 

(l-a)(l-£) = (l-e) (60) 

in which 

a = J2 a i > P = S— (61) 

i=1 i=1 P 
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a is the total probability of injection of particles (note that a should be less than 
one ) and (3 is the average probability of extraction of particles. In the special case of 
1-species (60) reduces to 

(l-a 1 )(l-/? 1 ) = l-d 1 

Comparing this with the usual ASEP [33] in which the condition for one dimensional 
representation reads to be (1 — cc)(l — f3) = 1—p ( p is the hopping probability), make us to 
take e as the average probability of hopping i.e. e = Yh=i ^- So a natural choice for diS 
would be to take them ^. In one dimensional representation, the hopping probabilities 
are restricted to 



a < vi < V2 < t>3.... < v p < 1 (62) 

Within one dimensional representation, the stationary state is uncorrelated and is 
given by \P S >= \p >® L where 



p >= - 

c 



V 

v 2 



C = S + -(V 1 + V 2 ...V P ) = S + -V 
p p 



(63) 



\ V T> 

The density and current of i-type particles are all site independent and are respectively 
given by equations (57) and (56) 



p«_(a,i) = 



2i 
p 



J<-(a, i) 



(64) 



e (I_l) + 2 u ^°> e (I_i) + 2 

One can define total density and the total current by summing over all kind of species 
and finds 



v 

p 



;l _i) + 2 



e(i-l) + 5 



(65) 



4.2 Infinite-species limit 
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At this stage we consider the limit p — > oo, and we assume that the hopping proba- 
bilities of particles are chosen from a continuous distribution P(v). Discrete quantities 
j;F(i) are transformed into f(v)P(v) and sums into integrals. Equations (64) and (65) 
take the form 



V(a,v)P(v) T , , vP(v) . . 

^•-■' Mh'^) J - ia ' v) = <i - dVpm (66) 

"P(ct) 

= e (i-i) + 1 >( a ) = ( (i-.|^H (67) 

where 

(1 — a)f P(v) r l v 

V(a,v) = } K ] and Via) = (1 - a) / — — P{v)dv (68) 

Although one has many choices for P(v), we first take the following [19]. It has the 
merit that T){ot) can be analytically evaluated. 

Pl(p)= (i^p (t, " a)m ' m -° (69) 

This is a normalized distribution that vanishes with some positive power in low- 
velocities and increases up to v = 1. The average hopping probability e is found to 
be 

rl ( m _|_ 1) 

e— vPAv)dv = -(1 — a) + a 

Ja w (m + 2) v 7 

expressing m in terms of e and a we have 

2e — a — 1 

»=^- (70) 
for m to be positive, (70) implies (e, a < 1) 

2e-a-l>0 (71) 

We first study the current-density relationship for a fixed hopping probability, e. In 
order to do this, we evaluate T>(a) with (68) and replace m from (70) 

J - (a ' e)= e(I-l) + f^ (72) 

/ 2e— 1— a 
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Figure 1: The current versus the density for different values of e in backward updating. 
Continuous lines refer to Pi(v) and dotted lines refer to P2(v). 

2e— 1— ae 

e) = — — j 2e - 1 ~° 1 - (73) 

' v ' ; e(± - 1) + 2 9 e ~ r V ; 

/ 2e— 1— a 

The above expressions gives the total current and total density in terms of two control 
parameters namely the total arrival probability a and the average hopping probabilty e. 
We now eliminate a between and numerically which then gives the current density 
diagram. This diagram is shown in figure (1) for two values of e. 

Remark: Total current and total density are in general functions of three 
control parameters e, a and m. Recalling that e is the average hopping probabilty , a is 
the total rate of injection and m determines the shape of hopping distribution function. 



15 



0.70 



0.75 



0.80 



0.85 



0.90 



0.95 



1.00 



P 



Figure 2: The current versus the density for different values of a in backward updating. 
Continuous lines refer to P\{v) and filled squares refer to P2(v). 

Equation (70) implies that only two parameters are independent. There is a one-to- 
one correspondence between the two dimensional parameter space defined by the surface 
(70) and the current-density space. versus p^ in fig (1) corresponds to ntersection 
of planes e = constant , with the surface defined by (70). We can instead look at the 
intersection of a = constant planes with the surface and find the corresponding curves in 

— p<_ plane. This is done by eliminating e between equations (72) and (73). Figure 
(2) shows these diagrams for some values of a. 

Finally we consider the curves of constant m in — p<_ plane. To obtain these 
curves, one should write and p^ in terms of a and m as follows 



a(a + m + l) 



(74) 



(a + m + l)(l-a) + a(m + 2)(1 + ^) 
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Figure 3: The current versus the density for different values of m in backward updating. 
Continuous lines refer to P\(v) and dotted lines refer to P2(v). 



a(m + 2)(l + ^) 

P ^ (a ' m) ~ (a + m+l)(l-a) + a(ro + 2)(l + £) (?5) 

Eliminating a between p<_(a,m) and </<_(«, m) would give us the current-density 
diagrams for a fixed value of m. Figure (3) shows these diagrams for some values of m. 
As can be seen, the current does not vanish at p<_ = 1. This can be explained by noticing 
that although at p<_ = 1, the chain is completely filled, still we have current via exchange 
processes. At = 1, the more decreasing m, the more approaches to zero. 

Using (72) and (74), we can also look at the behaviour of current itself as a function 
of control parameters. In figures (4) and (5), we show the dependence of on a ,e for 
some fixed values of e and a. Note that for each a, there is a lower limit of e which can 
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J 




a 



Figure 4: The current versus the arrival probability of particles for different values of e 
in backward updating. Continuous lines refer to P\{v) and dotted lines refer to Pi{v). 

be obtained through equation (70). 
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Figure 5: The current versus the total probability of hopping for different values of a 
in backward updating. Continuous lines refer to P\(v) and dotted lines refer to P2(v). 
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Our second choice of velocity distribution function is the following 



P2(v)= {m ^^ 2 2 \ v- a r(l-v) m>0 (76) 

It vanishes at v — a , v — 1 and has a maximum at v max = ^*rf ■ If m increases, v max 
approaches to one and if m decreases to zero , it approaches to a. Inserting p2(v) into 
(39) we arrive at 

3e - 2a - 1 

m = (77) 
using (67), (68) and (77), we express and p^ in terms of e, a and a, m 

ea(2a + 1 — 3e) 

6) = e(l - a) (2a - 3e + 1) + a(2ae - 3e + 1) ( ' 

a(2ae + 1 — 3e) , , 

P ^ (a ' 6j ~ e(l - a) (2a - 3e + 1) + a(2ae - 3e + 1) { ' 

a(2a + m + 1) 

J^(a,m) = - : ; --7T- r (80) 

V ; (l-a)(2a + m + l) + a(m + 3)(^ + l) V ; 

a(^ + 3)(— + 1) , x 

p( a , m) = i ' Km - J — -j- 81 

' V ; (l-a)(2a + m+l) + a(m + 3)(^ + l) V ; 

We now eliminate a between J«_(a, e) and p<_(a,e) which leads to current-density 
diagrams for fixed values of e. Dotted lines in fig (1) shows these diagrams for the same 
values of e . 

Similar to Pi(v), we can consider the current-density diagrams corresponding to constant 
a and m. These diagrams are shown by dotted lines in figures (2) and (3) respectively. 
Dependence of J<_ on a and e for P2{v) are also shown in figures (4) and (5) by dotted 
lines. Note that in figure (5), the curves obtained from P\{v) asymptotically approach 
to those of P2,{v). 

Here, we would like to reveal a feature of the infinite species limit which is somehow 
reminiscent of Bose-Einestein condensation [19]. Equation (68) implies that the density 
of particles with speed v is proportional to . Taking (69,76) for P(v) we have 

p(v) ~ v(v - a)" 1 ' 1 (82) 
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Recalling that a is the minimum speed of particles, equation (82) shows two different 
kinds of behaviour depending on whether m > 1 or m < 1 . 

I) If m — 1 > then p(v) — > for v — > a 

which means that density of low speed particles is small, i.e. most of the particles 
move with rather high speed. 

II) If m — 1 < then p{v) — > oo for u — > a. 

In contrast to the case I, here the density of low-speed particles are large and most of 
the particles move with low speed , which can be interpreted as appearing of the traffic 
jammed phase. 



5 p-Spesies ASEP with Forward Updating 

5.1 Formulation 

As stated in the introduction and section (2), instead of right to left (backward) 
updating, one can change the direction of updating and starts from the first site of the 
chain (forward updating) and updates from the left to the right in the same manner of 
backward updating. Most of the steps are similar to backward updating and we only 
write the results. The transfer matrix takes the following form 

= R L T L ^ L ...T^U (83) 

All the matrices are the same as in (17,18,19). The MPS for the steady state is 
written as [23] 

\P, >^=« W\A® A® ... <g> A\V » (84) 

Taking A and A to satisfy the same algebra (21-23), makes \P S >^ to be a stationary 
state i.e. T^\P S >^= \P S >_>. Here at first site i = 1 a "defect" A is created, then 
transmited forward until it reaches the last site i = L where it disappears. Next we 
consider formulas for the currents and densities. Here the situation is quite different and 
the difference between forward and backward updating reveals itself. The definition of 
currents reads from (49-51) and is replaced with (83). The mean current of 2-type 
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particles through site k is found to be 

< W\C L \V > 

Where jW is the same as equation (53), and C — E + Z)f=i A- We again demand that 
E and _Dj satisfy equation (33,34) which in turn let us revisit equation (42-45) and thus 
we have 

J (i) = diC (86) 
C = C (87) 



Putting (86,87) in (85) yields 



< 



<W\C L \V> v ' 

Also one can write the mean density of z-type particles at site k 

nf < W\C^HD t - ^\V > 
k <W\C L \V> 



5.2 One dimensional representation and infinite number of species limit in 
forward updating 

Again scaling all ZVs by a ^ factor, we now take and E to be c-numbers. Similar 

to backward update, they are — and £ respectively and the equations (58-61) remain 
the same. In one dimensional representation, the densities and the currents of z-type 
particles are all site independent and are respectively given by 

/Mo^) = e(1 P _ X) P v , J^(u,i) = e( i_ P 1) + v (9°) 

Comparing the above equations with their counterparts in backward updating, we see 
that currents do not change but forward density undergoes the following modification 
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J_>(q;, i) — J<-(a, i) — J(a, i) , P->(ot, i) — p<-(ct, i) — J (a, i) (91) 

The above relations reveals the difference between forward and backward updating. 
Similar relation between backward and forward densities is seen in [23]. We again define 
the total density and current by summing over densities and currents of all kind of species 



J^(a) = J_(a) = _ , = p < _(a) - J(a) (92) 

v a ' P 

Now we take the limit of p — > oo. Adopting the same distribution functions -Pi(v) , 
P2{v) and using (92), one easily can obtain and p^ as functions of e , a and m , both 
for P\{v) and -p2(^)- Similar to the backward scheme, the corresponding current-density 
diagrams can be obtained by eliminating one of the control parameters. These diagrams 
are shown in figures (6) to (8). 
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Figure 6: The current versus the density for different values of e in forward updating. 
Continuous lines refer to P\(v) and dotted lines refer to P2(v). 
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Figure 7: The current versus the density for different values of a in forward updating. 
Continuous lines refer to P\(v) and filled squares refer to P2(v). 
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Figure 8: The current versus the density for different values of m in forward updating. 
Continuous lines refer to P\(v) and dotted lines refer to P2(v). 
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Remark: 

Surprisingly as can be seen in fig(7), when p^ goes to zero, the value of J_> does not 
vanish. This is in contrast to other's results and is an exclusive effect appearing only in 
forward updating which can be explained by noting that, when the lattice is completely 
empty, in first site a particle is injected with the probability a and according to the 
multiplicative nature of the transition matrix is transfered through the lattice, hence one 
has a non-zero current. 

In general, the value of J_> at p^ = is equal to a and this point refers to the point 
(m = oo, e = 1, a) in parameter space. 

We would like to end this section with some remarks on sub-parallel updating scheme. 
In fact as stated in section 1, there are few exact results in parallel updating. The root of 
this difficulty is the non-local nature of transfer matrix which in contrast to the ordered 
sequential updating, can not be written as a product of local transfer matrices. A simpler 
case is to consider a sub-parallel updating scheme. In this scheme, one proceeds with two 
half time-steps. In the first half, one updates the first site, last site and all pairs (Tj, r i+ i) 
with an even i (L is taken to be even). Then in the second half time-step, one updates 
all pairs (r^, r i+ i) with i odd. So the transfer matrix is 

T sp = T$T$ (93) 

with 

T$> = LiT 2i3 T 4 ,5 • • • Tl^l^Rl (94) 

T$ =T 1 , 2 T 3 ,4---r L _ 1 , L (95) 
Defining MPS for sub-parallel updating as follows 



\P S > sp =«W\A(g) A(g) A(g) ■ ■ ■ A(g) A\V » (96) 

It can be verified that T sp \P s > sp = \P S > sp provided that equations (21-23) are satis- 
fied. 

It is shown in [32] that sub-parallel and ordered sequential updating schemes are 
intimately related to each other. It is proved that in general the following correspondence 
exists 

. (i) . f < k odd , _x 

< Wfc >s P = < f . (97) 

{ < n k >^ k even 
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/ (i) CO ^ f < nftnP >^ k , I odd , nQ , 
<nl ) n\>> > sp = f. } . } ' (98) 

[ < n y k 'n\ ' >^ k , / even 

where k and / refer to the lattice sites and i and j refer to the state of the site. 
Using this general correspondence, we obtain the density profile of p-species ASEP under 
sub-parallel updating ( one dimensional representation ) 

< 4° >s P = 7T _ P n ~ v k = even (") 

<4 ) >s P = ( sjJ.v k = odd ( 10 °) 

6 p-Species ASEP with Ordered Updating on a Ring 

In this section we consider the p-species ASEP on a closed ring of N sites. We work in a 
canonical ensemble in which the number of each species (i) is fixed to be rrii and we take 
the total number of particles to be M i.e. Ya=i m % — M . 

The periodic system can be described by a one dimensional representation of the bulk 
algebra (24-28). In this case the bulk algebra reduces to the following equations 

(1 - v t )die = d t e (101) 
(\^ i )d j d i = djdi (102) 



The above equations yield 

dj = - e (l - Vj )dj (103) 

Here di and di, correspond to one dimensional representations of Di and Di (not to 
be confused with those introduced in (34)). Using (53) and (57) we obtain the following 
forms for the density and the current of i-type particles : 

= J - = iM + ~ [vA ? d > ~ di ^ djVj ^ (104) 

e ^ l 1 e j j 

Summing over i, we obtain the total current and density 
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p^ = ^r^V ' J - = «E^ ( 105 ) 

Defining the population averaged velocity < v > as follows 

< v >= Y 106 
and rescaling the d^s and e so that 

e + J2di = e + J2di = l (107) 

i i 

we arrive at 

J_ = <v>P^-P-) (108 ) 
1 - < v >p^ v ; 

which is the current-density relation of jo-species ASEP on a ring with backward 
updating. Comparing it with the usual ASEP on ring with backward updating in [23] , 
we see that they both have the same form. In p species model, < v > plays the role of 
hopping probability in usual ASEP. Fig (9) shows versus p<_ for different values of 

< v >. 

The maximum current occurs at 

p-« v » = > I (109 ) 

H v ; < v > - 2 y J 

We now consider the forward updating. Note that since we don't have particle-hole 
symmetry, the current-density relation in forward updating can not be obtain from the 
one in backward updating and should be considered seperately. In forward updating we 
have 

(i) = j( ! ) = J W (110) 

Using (101-103) and (107), after straightforward calculations, we arrive at 

= 1 t LJ!L — i-f (ill) 

• l—v 

where 



29 




Figure 9: The current versus the density for different values of < v > in backward 
updating. 
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If we now take p — 1 , < > will reduce to and (111) takes the following form 

J_ = (113) 

1 - Vip^ 

and the particle-hole symmetry is recovered [23] i.e. (113) is obtained from (108) by 
changing to 1 — 

Fig(10) shows J_> versus p_> for different values of < >. 

The maximum of J_> has moved to the left. This maximum occurs at 

pZ ax {< >) = < iri; >)'-!] ^ ^ ( 114 ) 

l—v 
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Figure 10: The current versus the density for different values of < > in backward 
updating. 

7 Comparison and Concluding Remarks 

Here we compare our results with those of [20] and specify the similarities and differ- 
ences between ordered and random sequential updating procedures. Both procedures are 
described by a similar quadratic algebra. In random update, time should be so rescaled 
such that the average hopping rate equals one. On the contrary in ordered updating 
the average hopping probability e remains as a free parameter. This is one of the main 
differences between two updating schemes. In both schemes, injection rate (probability) 
of an z-type particle (ctj) is proportional to its hopping rate (probability) t>j. 
When considering infinite species limit, one can investigate the characteristics of both 
schemes with a limited number of control parameters. As long as analytical calculations 
are concerned, these control parameters are a , e and m in ordered schemes and a, m and 
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A in random scheme where m, A determine the shape of distribution function [20]. One 
of the advantages of ordered scheme is the appearance of the more physical parameter e 
in control parameters, which is absent in random scheme. 

In this paper, we made a more complete investigation of the current- density and current 
diagrams for different regions of parameter space. 

We also evaluated the dependence of the current on the density for fixed values of a in 
random scheme. The corresponding diagram is quite similar to ours in figure (2) but the 
values of current and minimum allowed value of the density are different. 
All the result of this paper and [20] have been obtained in a restricted region of pa- 
rameters space (ai,(3i,Vi) where mean field approximation becomes exact. It would be a 
highly nontrivial task to investigate the physical properties of the hole regions of parame- 
ter space either by infinite dimensional representations or by the explicit use of quadratic 
algebra. Consideration of p-species ASEP under fully parallel updating is another inter- 
esting subject which is under study. 
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